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^ ■ ABSTRACT 

\ High-redshift submillimetre-bright galaxies identified by blank field surveys at millime- 

tre and submillimetre wavelengths appear in the region of the IRAC colour-colour di- 
agrams previously identified as the domain of luminous active galactic nuclei (AGNs) . 
Our analysis using a set of empirical and theoretical dusty starburst spectral energy 
distribution (SED) models shows that power-law continuum sources associated with 
hot dust heated by young (< 100 Myr old), extreme starbursts at z > 2 also occupy 
the same general area as AGNs in the IRAC colour-colour plots. A detailed com- 
■ parison of the IRAC colours and SEDs demonstrates that the two populations are 

distinct from each other, with submillimetre-bright galaxies having a systematically 
flatter IRAC spectrum (> 1 mag bluer in the observed [4.5]-[8.0] colour). Only about 
20% of the objects overlap in the colour-colour plots, and this low fraction suggests 
\ that submillimetre galaxies powered by a dust-obscured AGN are not common. The 

red IR colours of the submillimetre galaxies are distinct from those of the ubiquitous 
foreground IRAC sources, and we propose a set of IR colour selection criteria for iden- 
, tifying SMG counterparts that can be used even in the absence of radio or Spitzcr 

^ ■ MIPS 24 urn data. 

' 
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1 INTRODUCTION 

One of the most exciting developments of the past decade 
has been the resolution of the cosmic far-infrared back- 
ground into discrete sources, providing a first glimpse of 
the rapid build-up of massive galaxies in the early uni- 
verse long predicted by theory. Deep, wide blank-field sur- 
veys at millimetre (mm) and submillimetre ( s ubmm) wave- 
lengths dSmail et al.lll997l; iBarger et al.lll998l;lHughes ; et all 



2006 
2005: 



Eales et al.l 1999; IScott et al.l 

Greve et al1l2004 [Laurent et al.l f2005: Mo rtier et all 



' 20021: IWang et alj|2004, 



iBertoldi et~ai1 120071 : IScott et all 120081 ) have shown 
that ultraluminous infrared galaxies (ULIRGs) at z > 1 
contribute significantly to the observed far-IR background. 
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Multi-wavelength follow-up studies of these so-called sub- 
millimetre galaxies (SMGs) suggest that they are massive, 
young galaxies seen during the period of rapid stellar mass 
build-up, with very high specific sta r formation rates at 
2 > 1 (see review bv lBlain et al.| [2002h 

A major obstacle in understanding the nature of this 
luminous dusty galaxy population is the limited angular 
resolution (> 10") of current instrumentation, which pre- 
vents unambiguous identification of their counterparts at 
other wavelengths. Deep interferometric radio imaging and 
S'pitzer 24 u.m MIPS imaging are shown to be effective 
for identifying a subset of SMGs at z < 3 (|lvison et al.l 
l2002t IChapman et al.l [20051 : iPope et alj|2006h . and our un- 
derstanding of their nature and evolution is based almost 
entirely on some 100 SMGs identified this way. Recent 
890 u.m continuum imaging of the seven brightest AzTEC 
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Wilson et ail20 08a) 1100 urn sources in the COSMOS field 



Scoville et al 



2007) using the Smith sonian Submillimeter 



Array (SMA) bv lYounger et al.l l|2007l ) has shown that there 
may be a substantial population of higher redshift (z > 3) 
SMGs that are extremely faint or undetected at radio and/or 
24 u.m MIPS bands. This discovery accentuates the need for 
a new method to identify and investigate the higher redshift 
SMGs in order to quantify their contribution to the cosmic 
energy budget and star formation history, and to map their 
evolution over time. 

One common feature among all seven COSMOS 
AzTEC/SMA so urces and t he fo u r other SMGs det ected 
with the SMA b vllono et all (|2006l l. lWang et all l|2007l ). and 
lYounger et al.l 1 20081 ') is that they are all detected in the 
Spitzer IRAC bands at > 1/xJy level, raising the exciting 
possibility that deep IRAC imaging may provide a powerful 
new tool for identifying and obtaining deeper understand- 
ing of the SMG phenomenon. The IRAC bands cover at 
z ~ 3 the rest-frame optical to near-IR portion of the spec- 
tral energy distribution (SED) and offer direct insight to the 
properties of their stellar component. A detailed analysis of 
the multiwavelength properties of the AzTEC/SMA sources 
is presented elsewhere (Yun et al., in prep.). In this paper, 
we examine the rest-frame optical/near-IR properties of a 
large sample of well-studied SMGs and the utility of the 
IRAC colour-colour plots for identifying and investigating 
the nature of the SMG population in general. 



2 Spitzer IRAC COLOURS OF 

SUBMILLIMETRE GALAXIES AND AGNS 

2.1 Red IRAC colours of SMGs 

A colour-colour diagram is a powerful tool for analyzing 
SEDs sampled coarsely with just a few broad-band mea- 
surements. The four Spitzer IRAC bands (3.6 u.m, 4.5 u.m, 
5.8 u.m, and 8.0 u.m) probe the portion of the galaxy SED 
that includes photospheric emission from cool stars (at 
z — — 4) and power-law continuum from hot dust sur- 
rounding young stars and/or AGN. Polycyclic aromatic hy- 
drocarbon (PAH) and other spectral features are also impor- 
tant for the z ~ galaxies. The 1.6 u.m stellar photospheric 
PAH feature is prominent in stellar systems older than 10 
Myr, and galaxies with substantial cool stellar populations 
appe ar blue in an IRAC colour-c o lour diagram a s a result 
(see ISimpson fc EisenhardtJ 1 19991 : ISawicki 120021 ). Sources 
dominated by stellar photospheric emission form a densely 
concentrated cloud near (—0.4, —0.4) in the S5.8/S3.6 vs. 
5*8. 0/54. 5 colo ur-colour diagram for a sample of field galax- 
ies shown bv lLacv et al.l (|2004T ) . The z ~ late type galax- 
ies with varying amounts of PAH emission in the 8.0 u.m 
band form a distinct branch of galaxies with a constant 
S5.8/S3.6 ratio emerging from this cloud. Noting that 54 
quasars identified in th e Sloan Data Release 1 quasar survey 
l|Schneider et alJl2003h are associated with a second branch 
with red IR colours, Lacy et al. identified a wedge-shaped 
area in the colour-colour diagram as the characteristic re- 
gion populated by obscured and unobscured AGNs. 

To explore the importance of AGN activity among 
SMGs, we examine the IR properties of a sample of 47 well 
studied SMGs and compare them to AGNs using the AGN 



diagnostic IRAC colour-colour plot by lLacv et al.l l|2004l ) . 
Our SMG sample is constructed from the literature primar- 
ily for their secure identification using high angular res- 
olution interferometric imagi ng at radio wav e lengths and 
deep ffpitzer MIPS imaging dlvison et all 20051 ; iGreve et al.l 
l2005l ; iTacconi et ail 120061 ; iPope et al.ll2006h . A small sub- 
set of this sample has be en observed with the ffpitzer In- 
frared Spectrograph (IRS ; Menendez-Delmcstrc et al. 2007; 
IValiante et"al] 120071 ; iRigbv et al.l 120081 ; IPope et all l2008h . 
and those classified as "starburst" (strong PAH features) 
and "starburst+AGN" (PAH plus power-law continuum) 
are shown using different symbols in Figure [1] Of these 
20 submm-selected galaxies, none shows a pure power-law- 
dominated IRS spectrum characteristic of an AGN. We com- 
pare the IR colours of our SMG sample to those of the 19 
z ~ 1 AGNs identified by their p ower-law continuu m in the 
First Look Survey (FLS) field by lLacv et all (|2004T ). and to 
the IR colours of the foreground galaxy population, using 
the 4000 I RAC sources random ly selected from the COS- 
MOS field (jSanders et alj|2007h . The IRAC colours of these 
foreground sources are slightly offset from the FLS centroid 
(indicated in Figure [1] by a large circle) , probably because 
the COSMOS sources are on average fainter and at higher 
redshift. 

A surprising result is that more than 90% of the SMGs 
are located with in the region designated for AGNs by 
lLacv et al l l|2004l ). It is tempting to interpret this result as 
an indic ation that th e majo rity of SMGs host a luminous 
AGN, as IGreve et all (120071 ) concluded for their MAMBO- 
selected SMGs in the field surrounding the z — 3.8 ra- 
dio galaxy 4C 41.17. However, the majority (13 out of 20, 
or 65%) of SMGs observed with Spitzer IRS show spectra 
dominated by PAH features, characteristic of pure starburst 
systems. The remaining seven SMGs show starburst+AGN 
hybrid spectra, suggesting that the power-law AGN emis- 
sion is not the dominant contributor to the total IR lumi- 
no sity of these SMGs e ither. Four of the SMGs targeted 
by IValiante et all l|2007l ) were undetected by their Spitzer 
IRS observations. Deep Chandra X-ray data in the GOODS- 
North field suggest that only a sm all fraction of SMGs ar e 
detected in the hard X-ray band (|Pope et al.ll2006l . 120081 ) . 
further supporting the notion that energetically-dominant, 
luminous AGNs are n ot common among the SMGs (also see 
lAlexander et al.ll2003l ). 

Another commonly used diagnostic for AGN activity is 
the Spitzer/IRAC [5.8]-[8.0 ] vs. [3.61— [4 . 51 colo ur-colour di- 
agram, first introduced by IStern et al.l l|2005h . In addition 
to red colours, Stern et al. utilize empirical colour tracks 
of star- forming galaxies (e.g., M82) to further differentiate 
starbursts from AGNs. We reproduce the Stern et al. plot in 
Figure [5] using the same samples as used in Figure [T] About 
1/3 of the SMGs (16/47), including 5 out of 13 "starburst" 
IRS spectrum sources, are now found outside the AGN re- 
gion outlined by Stern et al. However, a large fraction of the 
SMGs (31/47), including many with "starburst" IRS spec- 
trum, are within the AGN region. 



2.2 Analysis using empirical and theoretical SED 
models 

A natural explanation for the observed red IR colours 
of SMGs emerges when empirical and theoretical SED 
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Figure 1. Spitzer/IRAC S5.S/S3.6 vs. Sg. 0/S4, 5 colour— colour diagram for SMGs and IR QSOs. (a) Empty squares are SMGs se- 
curely identified by radio and CO interferometric imagin g. Filled squares and pentagons are SMGs identified as "starburst " and "star- 
burst +AGN" by Spitzer IRS spectra, respectively dMenendez-Delmestre et aljbool lValiante et alj|2007l; iRiebv et al.ll2008l ; iPope et all 
2008). IR QSOs with power-law spectrum identified in the FLS field by Lacy et al. and |MartjnezJ5ansj ereetH!dl I ||2008|'| arc shown as 
stars. Small dots represent 4000 random field galaxies selected in the COSMOS field llSanders et alj|2007l ), and the large circle cen tered 
near (—0.4,-0.4) represents the centroid of the IRAC sources in the First Look Survey (FLS) field identified bv lLacv et all j2004!). The 
long-dashed line outlines the region for AGNs proposed by Lacy et al. The solid line shows the extended region we propose for the 
identification of the SMG counterparts (Eq. 1). An error bar corresponding to a typical 10% uncertainty in the flux density measure- 
ment is shown on the right side, (b) The same Spitzer/IRAC colour-colour diagram for IRAC sources found within a 6" radius of the 
nine securely identified mm/submm sources. Filled circles are the nine SMGs identified by direct submm interferometric measurements 
obtained using the SMA while the empty circles represent the foreground/interloper I RAC sources, (c) Redsh ift evolution colour-colour 
tracks for three different starburst ages based on theoretical starburst SED models by lEfstathiou et alj {2000). Filled and empty squares 
along the redshift evolution colour tracks mark the redshifts of z = 0.5, f, 2, 3, 4, and 5. The thick solid line represents power-law 
spectrum sources with a=0. 3-1.0 (S v oc v~ a ). (d) The effects of extinction are demonstrated by the 72 Myr old starburst colour— colour 
model tracks with total extinction of Ay = 50 and 200. 



models of dusty starbursts are examined in the con- 
text of these IRAC colour-colour plots. The extreme 
IR luminosity of the mm/submm detected sources 
(Lir ^ 1O 12_13 L0) indicates that a dust-obscured ex- 
treme starburst or a luminous AGN dominates their en- 
ergy output, and the entire observed SE D can be de- 
scrib ed by a relatively simple model (see lYun fc Carillil 
2002). Theoretically motivated radiative tra nsfer mod- 
els for dust-enshrouded sta r bursts and AGNs llSilva et al 
19981; lEfstathiou et al.ll2000l; ISiebenmorgen fc KriigellbOO 



Chakrabarti fc Whitnevll2008l ) are successful in reproducing 
the observed UV-to-radio SEDs, including those of extreme 
objects such as SMGs and "hyperluminous" IR galaxies 



not commonly found in the local universe (see Farrah et alj 
2003; lEfstathiou fc Rowan- RobinsodlioollVega et al.ll2008l ; 



Groves et al.ll2008l ). A large number of free parameters (e.g., 



source geometry, initial mass function, metallicity) and some 
degeneracy among them limit the utility of these SED mod- 
els for unique quantitative interpretations of the photomet- 
ric data. Nevertheless, they are highly valuable tools for 
evaluating the effects of model parameters such as star- 
burst age and extin ction. For this wo r k, we adopt the model 
SEDs computed by lEfstathiou et alj (|2000T > for their ease of 
use and because they have already been widely tested (e.g. 
lEfstathiou fc Rowan-Robinsonll2003l ; iFarrah et al]|2003l ). 

The Efstathiou et al. SED models are shown in Fig- 
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Figure 2. Spitzer/TRAC [5.8]— [8.0] vs. [3.6]-[4.5] colour-colour diagram for submillimetre galaxies and IR QSOs in AB magnitudes, 
(a ) All symbo l s and models plotted are identical to those in Figure [T] The thin long-dashed line outlined the AGN region proposed 
by I Stern et al. | 11200511 . The solid line shows the extended region we propose for the identification of the SMG counterparts (Eq 2). (a) 
The same Spitzer/TRAC colour— colour diagram for IRAC sources found within a 6" radius of the nine securely identified mm/submm 
sources (filled circle). Empty circles are foreground/interloper IRAC sources found within the same search radius as discussed in the text. 
The thick-solid line represents po wer-law spectrum sourc es with q=0. 3-1.0 (S„ oc v~ a ) in both panels, (c) & (d) The same theoretical 
colour— colour evolution tracks by Efstathiou et al. (2000) shown in Figure[T]are plotted for comparison. 



ure |3] along with t he observed SEDs of the z = 2.49 SMG 
J123707. 7+621411 (|Pope et allboOrjl and the z = 1.74 QSO 
SSTXFL J172253.9+582955 (|Lacv et al.ll2004r ) in order to il- 
lustrate the IRAC colour evolution with starburst age. The 
shape of the model SED changes quickly during the first 100 
Myr, primarily driven by the rapid evolution of the young 
stellar population. As the stellar population ages, the com- 
bined effect of decreasing radiation intensity and increasing 
photospheric emission ("1.6 [im bump") makes the SEDs 
flatten (become bluer) in the IRAC bands. These are generic 
features of all SED models, only dependent on the details of 
the input stellar population synthesis model. Even though 
the build-up of cool stars and the 1.6 \im photospheric com- 
ponent become evident as early as ~30 Myr after the initial 
starburst, the overall colour of starburst systems remains 
red even at 64 Myr after the burst. In comparison, the red 
continuum of the IR AGN increases monotonically across 
the IRAC bands into the mid-IR (MIPS 24 urn) band, with 
a steeper slope than that of the SMG and most of the the- 
oretical starburst SEDs. The comparison shown in Figure [3] 



nicely demonstrates the clear difference in the origin of their 
near-IR emission and the spectral slope between IR AGNs 
and SMGs, despite the broad similarity in their red IRAC 
colours. 

Three colour tracks for a single starburst population 
SED model with different ages are shown in Figures []J & 
Oifc. They cover the full range of IRAC colours associated with 
both AGNs and SMGs, with the older starburst tracks show- 
ing successively bluer IRAC colours. The majority of SMGs 
appear scattered about the model SED colour tracks, con- 
sistent with their Spitzer IRS spectra being characteristic of 
starburst-dominated systems. For a given model SED, the 
IRAC colour becomes monotonically redder with increasing 
redshift at z > 1, and most SMGs have colours consistent 
with model SEDs redshifted to z = 1 — 5. The colour de- 
pendencies on starburst age and redshift are nearly parallel 
to each other, leading to some degeneracy between the two 
quantities. Nevertheless, the observed red IRAC colours of 
SMGs can arise only if SMGs are at high redshift (z > 1). 

The effects of extinction (column density along the line 
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Figure 3. A Comparison of the Efstathiou ot al. (2000) model 
SEDs with the observed S EDs of the z = 1.74 QSO 
SSTXFL J172253. 9+582955 l lLacv et alj 120041) and the z = 
2.49 SMG SMM J123707.7+621411 (GN19; iPope et al.|[200^) . 
The model SEDs are normalized to roughly match the mea- 
sured flux densities of GN19 in the IRAC bands. The SED of 
SSTXFL J172253.9+582955 is redshifted to z = 2.49 to match 
that of GN19 for an easier comparison. 



of sight) are not as important in determining the observed 
IRAC colours of SMGs, unlike the case at shorter, optical 
wavelengths. The model colour tracks for two different vi- 
sual extinctions (Ay = 50 and 200) shown in Figures [TJd 
& [2jl track each other closely at z < 2, despite the large 
opacity difference between the two models. Optical depth 
is greatly reduced at these long wavelengths (Ax /Ay = 
20 ~ 40 for the IRAC ba nds; see llndebetouw et all [20051 ; 
iRoman-Zuniga et al.l I2007T ), and the weak dependence on 
extinction Ay can be naturally understood. As redshift in- 
creases, the IRAC bands begin to probe the near-IR to op- 
tical bands, and an increasing dependence on extinction is 
expected. Indeed the IRAC colours of the two model SEDs 
diverge at 2 > 2, with the higher extinction Ay = 200 
model predicting redder IRAC colours as expected (see 
Fig. (TJI) , and some degeneracy between redshift and extinc- 
tion is also identified. Overall, extinction still plays a mi- 
nor role compared with starburst age, and we identify star- 
burst age and redshift as the dominant physical parameters 
that affect the observed IRAC colours. Accounting for SMGs 
with the reddest observed colours (e.g., log(S5.s/S , 3.6)>0.3, 
log(S , g.o/5 , 4.5)>0.3) requires young (< 30 Myr old) star- 
bursts at high redshifts (z > 3) or a power-law AGN domi- 
nating the rest-frame near-IR SED. 

An AGN-like warm IR colour is a generic feature for 
all stellar systems whose near-IR luminosity is dominated 
by a dust-obscured young stellar population. In analyz- 
ing the "mid-IR excess" of galaxies selec ted using t he In - 
frared Astronomy Satellite (IRAS) data, lYun et al.l (|200ll ) 
noted that the youngest dusty starbursts can display warm 
mid-IR colours, exce eding the classic Sey fert division at log 
(SWWSW) > 0.18 |de Griip et alj|l985l ). While exploring 



the presence of power-la w AGN candidates i n the Chandra 
Deep Field North region, iDonlev et ail (|2007h also noted the 
incursions by their empirical ULIRG colour tracks (derived 
from the observed SEDs of Arp 220, IRAS 17208-0014, and 
Mrk 273) into the AGN region in their IRAC colour-colour 
diagram. The incur sion of M82 - like o bjects into the AGN 
boundary has lead Stern et al.l 112 005) to refine the selec- 
tion boundary, and lBarmbv et al.l l|2006l ) question the com- 
pleteness and reliability of AGN identification using IRAC 
colours for a similar reason. Using theoretical and empir- 
ical starburst SED models, we demonstrate that young, 
dusty starbursts exhibit red IRAC colours, and red IRAC 
colour is not unique to power-law AGNs. The popular AGN 
identification met hods using red IRAC co lours, such as by 
lLacv et al.l (120041 ) and IStern et alj l|2005l ). should be used 
with caution and a clear understanding of this important 
caveat. 



2.3 Systematic colour difference between SMGs 
and AGNs 

Both SMGs and AGNs appear within the broadly defined 
red IRAC colour regions by Lacy et al. and Stern et al., 
but SMGs and AGNs are distinguished by clear system- 
atic difference in their mean IRAC colours. The SMGs as 
a group are systematically bluer than the AGNs, with only 
about 20% of SMGs appearing mixed in among the power- 
law AGNs in both diagnostic colour-colour plots. All AGNs 
show log(S'8.o/<S , 4.5) ^ 0.2 in Figure [1] while the overwhelm- 
ing majority (> 80%) of the SMGs have a smaller flux ratio. 
Similarly, all AGNs have [5.8]-[8.0] colours redder than +0.2 
in Figure [2] The most significant colour difference is between 
the 4.5 [im and 8.0 u.m bands with an average difference of 
([4.5] - [8.0]) > 1 mag. The [3.6]-[4.5] colour difference is 
the smallest, with ([3.6] - [4.5]) ~ 0.25 mag. 

The comparison of the full IRAC band SEDs shown 
in Figure [3] offers a visually compelling demonstration of a 
steeper spectral slope for the AGNs. The average IRAC SED 
of the AGN sample is in excellent ag reement with the em - 
pirical type-2 QSO SED template bv lPolletta et al] (|2007f ). 
The template SED for a type-1 QSO is nearly identical (not 
shown) to type-2 QSO at these wavelengths, and the steep 
spectral slope is characteristic of all obscured and unob- 
scured AGNs. In comparison, the average IRAC SED of the 
SMG sample is significantly flatter, with only a few out- 
liers with an AGN-like steeper spectral slope. The ULIRG 
Arp 220 SED redshift to z — 2 offers a far better match 
to the SMGs, lending further support for the dust obscured 
starburst interpretation. Furthermore, the flatter spectral 
slope (bluer IRAC colours) suggests that SMGs in general 
do not host a dust-obscured, energetically dominant AGN in 
most cases. We note that rigorous diagnostics of an AGN in- 
clude the detection of high-ionization, high-excitation emis- 
sion lines and copious amounts of hard X-ray emission, and 
an IR power-law spectrum is only indirect evidence for the 
presence of an AGN, and does not directly measure the AGN 
accretion power. 

The distinct difference in the observed SEDs between 
AGNs and SMGs also offers an important constraint on 
the possible SMG-QSO evolutionary scenario. The ap- 
parent correlation observed between the black hole mass 
(Mbh) and the host spheroid mass (velocity dispersion 
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Figure 4. The IRAC band continuum spectral energy distribu- 
tions for the AGNs and SMGs shown in Figures [JJ & [2] All SEDs 
are normalized at 4.5 um band, and the AGNs and SMGs are 
displaced vertically by an arbitrary amount for ease of compari- 
son. The 2 = 2 SEDs of typ e-2 QSO and Arp 220 (M82 SED is 
nearly identical) compiled bv lPolletta et all l|2007l ) are shown for 
comparison, with an arbitrary offset in the flux scaling. 



a: see iMagorrian et all Il998l : iFerrarese fc Merrittl l200d : 
iGebhardt et al l l2000l ) has raised considerable interests in 
the process that leads to the build-up of the stellar mass in 
galaxies and the growth of the central super-massive black 
hole. The evolutionary scenario that a massive nuclear star- 
burst associated wi th an ultralum i nous i nfrared galaxy leads 
to a QSO phase (jSanders et all 1 19881 ; H orman fc Scovilld 
1988) offers an attractive explanation for the Mbh — o re- 
lation if the SMG phase corresponds to the period of rapid 
stellar mass build-up. However, only the minor overlap be- 
tween the AGNs and the SMGs we find here suggests that 
the duration of the transition period should be much shorter 
than the SMG or the IR AGN phase in such an evolutionary 
scenario. 



3 IDENTIFICATION OF OPTIC AL/IR 

COUNTERPARTS TO SUBMILLIMETRE 
SOURCES 

As discussed above, every SMG identified with high angu- 
lar resolution interferometric imaging using the SMA has a 
clear IRAC counterpart, including those undetected in the 
radio and in the MIPS 24 u.m bands. S'pitzer IRAC data 
may therefore offer a powerful new method for identifica- 
tion of SMGs that is otherwise extremely difficult because 
of their faintness at optical and near-IR wavelengths. The 
uncertainty in the mm-wave positions of SMGs is typically 



5-10", and the high IRAC source densit y (~ 60 arcmin -2 
at the 1.4 /^Jy level in the 3.6 urn band; [Fazio et al1l2004r ) 
is too high to enable unique identifi cation of SMG co unter- 
parts in general (see discussions by IPope et alj [2006) . The 
new knowledge obtained from Figures [1] & [2] that SMGs as 
a population have red IR colours, similar to AGNs and dis- 
tinct from the foreground field population, offers the exciting 
possibility of identifying or significantly narrowing down the 
counterpart candidates using the IRAC data alone. 

The distribution of SMG IRAC colours and the star- 
burst model SED colour tracks in Figures [1] & [2] extends 
slightly beyond t he A GN selection boun daries proposed by 
lLacv et al.il|2004h and lStern et al.l l|2005h . We propose a new 
set of SMG counterpart candidate selection criteria for the 
S5.S/S3.6 vs. SS.0/S4.5 colour-colour diagram, expanded to 
include all SMGs and the model SED tracks as shown with 
long-dashed lines in Figure[T] The new proposed criteria are: 

log(S a .o/SiM) > -0.3 A log(S 5 .a/S 3 .6) > -0.3 

A log(Ss.o/S4. 5 ) < log(S 5 . 8 /S 3 . 6 ) + 014 

where A is the logical AND operator. Equivalent criteria 
for the Stern et al. colour-colour diagram (Fig. [2} in AB 
magnitudes are: 



([5.8] - [8.0]) 



> 
A 
A 



-0.4 
([3.6] 
([3.6] 



[4.5]) > 0.036 x ([5.8] - [8.0]) - 0.318 
[4.5]) > 2.5 x ([5.8] - [8.0]) - 2.5 (2) 



The relative merits of these selection criteria are discussed 
further below. 

To evaluate the utility of these new IRAC colour selec- 
tion criteria for the SMG counterpart candidate selection, we 
examine the colours of the IRAC sources in the fields sur- 
rounding the nine SMGs securely identified through submm 
interferometric imag i ng jlono et al.ll2006l ; IWang et al]|2007l ; 
I Younger et ai1l2007l . I2OO8I ). The 20" x 20" IRAC 3.6 um 
images (not shown) centered on the original submillimetre 
source coordinates contain up to a dozen IRAC sources each, 
demonstrating the difficulty of using the IRAC source cat- 
alog alone for the SMG counterpart identification. There 
are a total of 20 IRAC sources located within a 6" radius 
(typical positional accuracy) of the nine nominal SMG posi- 
tions; their colours are shown in Figures[Tp &[2]3. Colours of 
the nine SMGs, shown as filled circles, fall within our pro- 
posed SMG colour boundaries in both plots. These SMGs 
were not used in defining the IRAC colour selection crite- 
ria and they thus represent an independent and success- 
ful test of our method. The remaining IRAC sources have 
colours consistent with the background sources, and about 
1/2 of them fall within the new SMG candidate identifica- 
tion boundaries. The new colour selection technique does 
not completely resolve the confusion problem. However, the 
situation is now greatly improved, as the candidate counter- 
parts are reduced to a single unique candidate or two po- 
tential candidates, making expensive follow-up observations 
more efficient. Other identification methods such as deep 
radio or MIPS 24 \im imaging do not greatly improve the 
situation in these test cases since many of these SMGs are 
undetected at other wavelengths. Three of the securely iden- 
tified sources have IRAC colours typical of z ~ 2 starburst 
systems while the remaining six cluster near the power-law 
AGNs (see Fig. [TJj &[2p). In addition to being undetected 
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in the radio, these six SMGs with AGN-like colour are un- 
detected in the S'pitzer MIPS 24 [im bands, ruling out the 
presence of a power-law AGN component that extends into 
the mid-IR bands, as seen for the QSO spectrum shown 
in Figure El Therefore these SMGs with very red IRAC 
colours may be the youngest, highly-obscured starburst sys- 
te ms at very high reds hifts (z > 3; also see the discussions 
in lYounger et alj|2007| ). 

The IRAC colour selection method described here is 
qualitatively similar to the IRAC 8.0 u.m band selection 
method used to identify the SCUBA galaxies in the CUDSS 
14 hour field bv lAshbv etail (|2006l ). Noting the red colour 
of the SMGs and reduced foreground confusion compared 
with the shorter wavelength IRAC bands, Ashby et al. ar- 
gued that 8.0 |J.m selection offers a better means for iden- 
tifying SMG counterparts than near-IR or optical selection. 
Their 8.0 \im selection resulted in a different counterpart 
from the optical or if -band selection methods in ten out of 
17 cases, casting some doubt on the earlier identifications of 
SMGs. The main difference of the new IRAC colour selec- 
tion method outlined here is that we fully quantify the red 
colour of the SMG population using Eqs. 1 & 2 for easy im- 
plementation and calibrate them using a large, well-studied 
sample in conjunction with empirical and theoretical dusty 
starburst SEDs. 



4 DISCUSSION AND CONCLUDING 
REMARKS 

High-redshift submillimetre-bright galaxies identified by 
blank field surveys at millimetre and submillimetre wave- 
lengths appear in the region of the IRAC colour-colour dia- 
gr ams previous l y iden tified as the doma i n of lu minous AGNs 
by lLacv et"ai] |2004h and IStern et all (|2005r ). Rather than 
interpreting this as a sign that the majority of the SMGs 
are powered by a luminous AGN, we have shown using 
empirical and theoretically motivated dusty starburst SED 
models that their IR colours can be interpreted as those of 
power-law continuum associated with hot dust heated by 
young (< 100 Myr old), extreme starbursts at z > 2. These 
SMGs fall along the branch of galaxies extending from the 
blue photospheric peak toward the red power-law region in 
the IRAC colour-colour plot by lLacv et al.l l|2004l ). and our 
analysis suggests that this branch is a heterogeneous en- 
semble of power-law AGNs and dust obscured starbursts 
whose continuum near-IR luminosity exceeds that of the 
photospheric emission from cool stars. In fact, our analy- 
sis demonstrates that the popular red IRAC colour selec- 
tion methods for AGN identification, such as the criteria by 
Lacy et al. and Stern et al., should be used with caution 
because of the significant starburst contribution expected. 
While there is some overlap between SMGs and AGNs in 
these IRAC colour-colour plots, SMGs are systematically 
bluer ({[4.5] - [8.0]) > 1 mag), consistent with 30 to 70 
Myr old starbursts observed at redshifts between z ~ 1 and 
z ~ 5. 

Our examin a tion o f the model dusty starburst SEDs of 
lEfstathiou et al. I <|2000h shows that the main physical pa- 
rameters that determine the observed IRAC colours are 
starburst age and redshift, while extinction (column den- 
sity) plays a less important role. For the models examined 



in Figures \T\ and [2] starburst age and redshift are nearly 
degenerate. The use of I RAC colours as a redshift in dica- 
tor was f irst p roposed bv lSimpson fc EisenhardU l|l999h and 
ISawickil (|2002l ) . and various empirical photometric redshift 
relations have been proposed recently w ith estimated ac- 
curacies of 5z/(l + z) = 0.1 - 0.2 (see iPope etafl 120061 ; 
IWilson et al . 2008b). This photometric redshift method may 
not work well for galaxies with red power-law continuum 
and/or weak 1.6 u.m photospheric feature that are more 
common at z > 2. The degeneracy between the starburst 
age and redshift found with these theoretical model SEDs 
suggests that a careful calibration using a large sample of 
SMGs is needed in order to quantify fully the systematic 
uncertainty of this method. 

A detailed comparison of the IRAC colour-colour plots 
and SEDs shows that AGNs and SMGs are distinct from 
each other due to intrinsic differences in their energy source 
and dust distribution. SMGs as a group have a flatter SED 
(bluer by ([4.5] — [8.0]) > 1 mag) in comparison with AGNs. 
Only 20% of the objects overlap in the colour-colour plots 
shown in Figures Q] & [5] and this suggests that SMGs 
powered by an AGN are not commo n. In the context of 
the ULIRG-QSO evolutio nary scenario (|Sanders et al.|[l988l ; 
iNorman fc Scovillel I 1988h . the little overlap between the 
AGNs and the SMGs population may indicate that the tran- 
sition period is much shorter than the duration of the SMG 
or the IR AGN phase. 

The red IR colours of the SMGs are distinct from the 
colours of the majority of the foreground IRAC sources, 
which are both early and late type galaxies at z < 2. We 
show that colour selection criteria similar to those of AGNs 
proposed by Lacy et al. and Stern et al. can be used to pin- 
point the IRAC counterpart to an SMG uniquely or to pare 
down the candidates to just a few, improving the efficiency 
of expensive follow-up observations at other wavelengths. 
The lack of high-quality multi-wavelength data for a large 
sample of SMGs is currently the primary limiting factor that 
prevents a better understanding of the SMG phenomenon. 
The IRAC colour selection method discussed here appears to 
work even in cases when the radio and MIPS 24 \ixa counter- 
parts are too faint to be detected at the present sensitivity, 
and an unbiased investigation of the entire SMG population 
may become possible using this new identification method. 

There are some important advantages and disadvan- 
tages to using the (Ss.s/Sa.e vs. Ss. 0/S4.5) colour-colour plot 
(Figured]) versus the ([5.8]-[8.0] vs. [3.6]— [4.5]) colour-colour 
plot (Figure[2]) for SMG candidate identification. The former 
uses IR colours with a longer baseline in wavelength, making 
the colour measurements more robust and the model colour 
tracks better behaved, as shown in Figure [1] On the other 
hand, the instrumental sensitivity of the IRAC 3.6 u.m and 
4.5 [im bands are a factor of a few better than the 5.8 u.m 
and 8.0 nm bands, and the 5.8 [im and 8.0 u.m channels will 
be inoperable during the upcoming warm S'pitzer operation. 
The model colour tracks shown in Figure [2] are more irregu- 
lar because of shorter wavelength baselines in these colours, 
but the [3.6]-[4.5] colour is the best determined quantity 
among all IRAC colour combinations. 

A far simpler and potentially more robust alternative 
SMG candidate selection criterion is [3.6]-[4.5] > —0.2. This 
single colour selection method does only a slightly poorer 
job of rejecting foreground sources than the full colour sec- 
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tion criteria described by Eq. 2. Given the superior sen- 
sitivity of the IRAC 3.6 u.m and 4.5 |i.m bands and their 
availability during the warm Spitzer mission, this simpler 
colour selection criterion may be the more effective method 
for identifying SMG counterpart candidates in the longer 
term. Alternatively, colours of IRAC to optical or near-IR 
bands or to MIPS 24 [im bands have also been proposed pre- 
viously for identifying AGN activity (see iHuang et aUl2004l ; 
I Ashbv et ail 12006k IWebb et all 120061 ). and similar analysis 
may also prove fruitful for SMG candidate identification. 
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